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ABSTRACT 

We present the characterization of the star KOI 961, an M dwarf with transit signals indicative of 
three short-period exoplanets, originally discovered by the Kepler Mission. We proceed by comparing 
KOI 961 to Barnards Star, a nearby, well-characterized mid-M dwarf. By comparing colors, optical 
and near-infrared spectra, we find remarkable agreement between the two, implying similar effective 
temperatures and metallicities. Both are metal-poor compared to the Solar neighborhood, have low 
projected rotational velocity, high absolute radial velocity, large proper motion and no quiescent Ho 
emission-all of which is consistent with being old M dwarfs. We combine empirical measurements of 
Barnard's Star and expectations from evolutionary isochrones to estimate KOI 961's mass (0.13 ± 
0.05 Mq), radius (0.17 ± 0.04 R ) and luminosity (2.40 x i()- 30±0 - 3 L ). We calculate KOI 961's 
distance (38.7 ± 6.3 pc) and space motions, which, like Barnard's Star, are consistent with a high 
scale-height population in the Milky Way. We perform an independent multi-transit fit to the public 
Kepler light curve and significantly revise the transit parameters for the three planets. We calculate 
the false-positive probability for each planet-candidate, and find a less than 1 % chance that any one of 
the transiting signals is due to a background or hierarchical eclipsing binary, validating the planetary 
nature of the transits. The best-fitting radii for all three planets are less than 1 R®, with KOI 961.03 
being Mars-sized (Rp — 0.57 ± 0.18 Re), and they represent some of the smallest exoplanets detected 
to date. 

Subject headings: Stars: individual (Barnard's Star, KOI 961) — Stars: low-mass — Stars: funda- 
mental parameters — Stars: late-type — Planetary Systems 



1. INTRODUCTION 

In February of 2011 the Kepler Mission announced 997 
stars that show light curv es consistent with t he pres- 
ence of transiting planets (|Borucki et al.l 120111 ). Many 
of the Kepler Objects of Interest, or KOIs, host planet- 
candidates that are terrestrial-sized, and represent 
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some of the smallest exoplanets det ected around main- 
sequen ce stars if they are confirmed. Morton fc Johnson! 
(|2011[ ) found that many of the planets have low a priori 
probabilities of being false positives, and careful follow- 
up studies will validate the planetary nature of many of 
the transit signals. 

KOI 961 {Kepler ID: 8561063) is particularly interest- 
ing: an M dw arf star with large pro per motion (/i = 431 
± 8 mas/yr, iLepine &: Sharal 120051) and transit signals 
indicative of three exoplanets, all with orbital periods 
of less than 2 days. The inferred radii of the planet- 
candidates depends directly on the assumed radius of 
the host star. Unfortunately, the stellar param e ters o f 
KOI 961 are highly uncertain. IBorucki et al.l (|2011[ ) 
listed an effective temperature of 4188 K and a radius 
of 0.68 R© based on pho tometric analysis fro m the Ke- 
io/er Input Catalog (KIC, iBatalha et al. 1120101 ). However, 
IBrown et alj (|2011[ ) warned that the stellar radii listed in 
the KIC have uncertainties of 30 to 50%, and that the 
radii for M dwarfs are particularly untrus tworthy. As- 
teroseismology studies of Kepler objects bv lVerner et al.l 
(|2011f ) also indicate biases in the surface gravity and ra- 
dius estimates in the KIC. 

Typically, M dwarf masses and radii are estimated 
by combining stellar lu minosity with reli able mass- 
luminosity relations (e.g. Deifo sse et al.l[2000| ) and mass- 
radiu s relations predicte d by stellar evolutionary models 
(e.g. iBaraffe et aLI 119981) . often with corrections to ac- 



2 



count for discrepancies between measured and predicted 
radii (e.g. iTorresI I2007D . Unfortunately, KOI 961 does 
not have a parallax measurement in the literature, which 
is necessary to estimat e its luminosity and hen ce mass. 

In a previous paper (jMuirhead et al.ll2011bl hereafter 
Paper 1), we acquired a if -band spectrum of KOI 961 as 
part of a survey of low-mass KOIs. The spectrum con- 
firmed that it is a dwarf and not a giant based on the 
shallowness of the CO absorption band, and we reported 
an effective temperature of 3200 ± 50 K and an over- 
all metallicity ([M/H]) of -0.33 ± 0. 12 using the spec- 
trosco pic indices and calibrations of iRojas-Ayala et al.l 
(|2011l hereafter RA11). KOI 961 was the coolest M 
dwarf in our survey, and it is likely the least massive 
KOI in the February 2011 Kepler data release. We in- 
terpolated the effective temperature and metallicity of 
KOI 961 ont o the 5-Gyr Solar-m etallicity evolutionary 
isochrones of lBaraffe et al.l |l998), providing a stellar ra- 
dius estimate of 0.19 R© and a mass estimate of 0.16 M Q . 
The stellar parameters derived for the low-mass KOIs in 
Paper 1 are spectroscopic, and are more accurate than 
the KIC values, but depend strongly on the accuracy of 
evolutionary isochrones for low-mass stars. 

In this paper we present measurements of the stel- 
lar parameters of KOI 961 by direct comparison to a 
well-studied low-mass star: Barnard's Star, or Gl 699. 
Barnard's Star is the second nearest main-sequence star 
system to the Sun at a distance of 1.824 ± 0.00 5 pc 
(jvan Leeuwenl l2007'). Discovered bv [Barnard! (|1916l ). his 
"Small Sta r with Large Prop er Motion" is classified as an 
M4 dwarf (jReid et al.lll995f) . and lLeggettl (j 19921 ) showed 
that Barnard's Star's kinematic and photometric prop- 
erties are consistent with an "old disk-halo" star. It has 
low projected rotatio nal velocity (V sm(i) < 2.5 km s -1 , 
iBrowning et al.ll201Q[ ). infreque nt Ha flaring with n o Ha 
emission during quiescence (e.g. lPaulson et al.ll2006fl and 
lies b elow the average main sequence ([Giampapa et al.l 
1996); all of which are consistent with an old, metal- 
poor M dwarf. Being nearby, Barnards Star is among the 
most studied and best characterized M dwarfs despite its 
extreme dissimilarity to the Sun. It has an empirically- 
measured s tellar radius usin g optical long-baseline inter- 
ferometry (lLane et al.ll200lD, a carefully measur ed bolo- 
metric luminosity (|Dawson fc De Robertiil2004|) , and an 
accurate Hipparcos parallax, which, when combined with 
photometry and the empiri c ally-d erived mass-luminosity 
relations of Delfos se et al.l |2000), provides an accurate 
measurement of its mass. 

RA11 report metallicities for Barnard's Star of [M/H] 
= -0.27 ± 0.12 and [Fe/H] = -0.39 ± 0.17, and an ef- 
fective temperature of 3266 ± 29 K. The same methods 
were used to measure similar values of KOI 961. This 
motivates a side-by-side comparison of the two objects, 
for the purpose of measuring the stellar parameters of 
KOI 961. In Section^ we perform a qualitative compar- 
ison of the two stars' spectra, finding them to be very 
similar. We also compare various spectroscopic indices 
and argue that KOI-961 is slightly metal-poor compared 
to Barnard's Star, yet shares its effective temperature, 
consistent with the K-band measurements. 

In Section [3] we perform a more quantitative assess- 
ment of KOI-961's stellar properties relative to Barnard's 
Star. We start with empirical measurements of the mass, 
radius and luminosity of Barnard's Star and appeal to 



stellar evolution models to estimate the relative differ- 
ence between the two stars' physical parameters. In this 
manner we find that KOI-961 is slightly smaller and less 
massive than Barnard's Star. We also estimate KOI- 
961 's distance and space motion and show they are con- 
sistent with a high scale-height, or "thick-disk," popula- 
tion in the Milky Way. 

In Section|4]we refit transit curves for the three planets 
using the public Kepler data, revise the transit param- 
eters, and revise the planet parameters, finding that all 
three are sub-Earth-sized. In Section [5] we calculate the 
false-positive probability for the transit signals based on 
a priori statistics of eclipsing binaries and seeing-limited 
images, historical images, and the high-resolution opti- 
cal spectrum of KOI 961, and we validate the planetary- 
nature of the transits. In Section|6]we discuss the masses 
and dynamical configuration of the three planets. In Sec- 
tion [7] we discuss the implications given that KOI 961 is 
metal-poor, and one of only a few mid-M stars surveyed 
by the Kepler Mission. 

2. COMPARISON OF OBSERVABLES 
2.1. Photometry 

Digital photometric measuremen ts of KOI 961 include 
griz photometry from the KI C dBatalha et al. I2010T) 
and 2MASS JHK photometry (jCutri et al.l 120031) . Dig- 
ital photometric measurements of Barnar d 's Sta r in- 
clude UVBRI photometry by iKoen et all (|2010t ) and 
2MASS JHK photometry. The 2MASS J- and F-band 
magni tudes of Barnard's St ar are near the saturation 
limit (Skrutsk ie et al.l 12006:. and produce colors that 
are in consistent with the spectral type from iReid et al.l 
(|1995| ) . With this in mind, we used the Jqtt- and Hqat - 
band measurements of Barnard's Star bv lLegg ett (1992), 
which we converted to 2M ASS J and H ma gnitudes us- 
ing the transformations of iCarpenterl ()200lD . 

We obtained Johnson B- and F-band and Cousins Rc- 
band measurements of KOI 961 with the Hereford Ari- 
zona Observatory (HAO) 11-inch telescope on UT 2011 
Novem ber 17, using standar d star fie l ds fro m iLandoltl 
(|1992f ) for calibration. Three ILandoltl (|1992t ) star fields 
were used containing a total of 45 total stars with accu- 
rate B- and F-band photometry, and 40 stars with ac- 
curate i?c-band photometry. We acquired 108 B- and 
y-band measurements of reference stars, and 94 Rq- 
band mea surem e nts to calibrate the measurements of 
KOI 961. I Gary! (|2007| ) describes the calibration proce- 
dures for HAO photometric measurements. The mea- 
surements and uncertainties are listed in Table [TJ along 
with available 2MASS J, H and Kg measurements. 

Figures Q] and [2] show J-Kg vs V-Kg and J-H vs. H- 
Kg for a sample of nearby M dwarfs, with Barnard's 
Star and KOI 961 indicated. The V-Kg color is sensitive 
to both effective temperatur e an d metallicity dLeggetti 
[1991 1 Johnson fc Appsl [2009h . In I Johnson et al.l (|2011l 
hereafter Paper 2) we identified the metallicity sensitivity 
of J-Kg. The J-H vs. H-Kg color-color regime has also 
been shown to distinguish metallicities and ages of M 
dwarfs, and w e include the population designations of 
iLeggettl (|1992D in Figured young disk (YD), young-old 
disk (Y/O), old disk (OD), old disk-halo (O/H) and halo 
(H). 

Barnard's Star and KOI 961 show remarkable agree- 
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TABLE 1 
Photometry of KOI 961 



Band 


Magnitude 


System 


B 


17.96 ± 0.14 


Johnson 1 


V 


16.124 ± 0.055 


Johnson 1 


R C 


15.055 ± 0.051 


Cousins 1 


J 


12.177 ± 0.021 


2MASS 2 


H 


11.685 ± 0.018 


2MASS 2 


K S 


11.465 ± 0.018 


2MASS 2 



"Obtained with the HAO 11-inch telescope. 
toutri et all 1(2003 ) : ISkFutskie et al.l (2006) 
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Fig. 1. — Color-color plot showing J-Kg vs V-Kg for stars 
within 20 pc with V — K > 3.0 (spectral ty pes later than K7) , 
trigonometric parallax uncertainty of < 5 % llvan Leeuw enl 120071 ; 
Ivan Altena et al.ll200lTl . V error of < 0.1 IIKoen et al.M20ia) and 
2MASS J and K s uncertainties of < 0.05 llCutri et alj 12003 ; 
Skrutskie et al. 2006). Note the difference in scales for the two 
axes. Included are Barnard's Star {black circle) and KOI 961 {red 
square), which are both at the lower edge of the distribution, ow- 
ing to their low metallicity. The similar V-Kg and J-Kq colors for 
the two stars indicate similar effective tempe rature and metallic- 
ity (e.g. Leggctt 1992; Johnson & Apps 2009), consistent with the 
A'-band measurements in Paper 1 and by RAIL 

ment in both color-color plots. This indicates that KOI 
961 and Barnard's Star have similar effective temper- 
atures and metallicities, and are metal-poor compared 
to the Solar neighborhood, consistent with the if-band 
spectroscopic measurements in Paper 1 and by RAll. 
KOI 961 has a slightly bluer V-Ks color (0.33 ± 0.04), 
which can be due to either higher effective temperature 
or lower metallicity. Medium-resolution optical spec- 
troscopy of the the stars indicates the latter, and we 
discuss this in the following section. 

2.2. Medium-Resolution Optical Spectroscopy 

We observed KOI 961 on UT 22 July 2011 with the 
Double Spectrograph on the Palomar Observatory 200- 
inch Hale Telescope: a two-chann el (Red/Blue), low - 
to-medium resolution spectrograph ()Oke fc Gunnl [l982 h 
We used settings that provided spectral coverage from 
5300 A to 10000 A in the red channel, with 10 A-wide 
resolution elements. At the time, the efficiency of the 
red channel limited the usable spectral coverage to 5300 
to 8000 A. The blue channel exposure containes signifi- 
cantly less signal, so we do not include those data. We 
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Fig. 2. — Color-color plot showing J-H vs. H-Ks for nearby 
M dwarfs, Barnard's Star {black circle) and K OI 961 {r e d squ are). 
As in Figure [T] we use J and H values from ILeggetli l|1992j ) for 
Barnard's Star out of concern for sat uration effects. We include 
the population designations of Leggctt (1992), which relate to the 
spread in metallicity: young disk (YD), young-old disk (Y/O), old 
disk (OD), old disk-halo (O/H) and halo (H). Both KOI 961 and 
Barnard's Star are in old and metal-poor populations. 

observed the white dwarf standard BD +284211 at the 
same airmass to calibrate the throughput and telluric 
absorption in the raw spectru m. We used a theoretical 
spectrum of BD +284211 from lStond (|1977D for the cal- 
ibration. We used available dome lamps for fiat-fielding 
and wavelength calibration, and reduced the data using 
the NOAO doslit package, written in IRAF. The final 
spectrum has a signal-to-noise ratio of roughly 350 at 

6 500 A. 

iDawson fc De Robertisl J2004) combined spectra of 
Barnard's Star taken over many wavelength regimes in 
order to accurately measured the star's bolometric lumi- 
nosity. We show a region of their combined spectrum 
and our spectrum of KOI 961 in Figure |3l The spec- 
tra are strikingly similar, aided by the fact that they 
have similar resolution. To quantitatively compare them, 
we calculated the spectroscopic indices of the Palomar- 
Michigan State U niversity Nearby - Star Spectroscopic 
Survey (PMSU e .g. IReid et aLlfilM IHawlev et al1H997t 
iGizis et aT1l2002|) . We used the band definitio ns and in- 
dex formula introduced bv IReid et al.l (fl995 h who also 
report the indices for Barnard's Star 

sp ectros c 

Barnard's Star from IReid et al.l (1 19951). o ur calculations 
using the IDawson fc De Robertisl (|2004D spectrum of 
Barnard's Star, and our calculations using our spectrum 
of KOI 961. The difference bet ween our ca l culate d in- 
dices for Barnard's Star and the IReid et all ()1995l ) val- 
ues is likely due to the spectra having different resolution. 
The standard deviation in the difference between our cal- 
culated indices of Barnard's Star and the PMSU indices 
is 0.05, which we used to estimate the uncertainty in the 
indices. 

KOI 961 has very similar CaH indices to Barnard's 
Star, but systematically higher TiO indices, indicating 
weaker TiO absorption. TiO absorption has long been 
the standard for spectral typing of M dwa rfs, and is 
sensi tive to stella r effective temperatu re (e.g. IReid et all 
1995). However, iLepine et al.1 (|2007f ) showed that TiO 



Table [2] lists the PSMU spectros copic indices for 
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Fig. 3. — Comparison of moderate-resolution optical spectra of 
KOI 961 and Barnard's Star. Red: Double Spectrograph spec- 
trum of KOI 961. B l ack: Spectrum of Barnard's Star from 
Dawson & De Robcrtis (2004). The spectra were renormalized, 
and the spectrum of Barnard's Star is offset for comparison. The 
stars have remarkably similar spectra indicating similar spectral 
type and metallicity. To quantitatively compare the stars we mea- 
sured the PMSU spectroscopic indices (see Tableland Figure [4j. 
Neither spectrum shows Ha e mission, but Barnard 's Star is known 
to have infrequent Ha flaring tPaulson et al.ir2006l) . 

is sensitive to both temperature and metallicity, and 
showed that CaH is a better diagnostic for effective tem- 
perature and spectral type. The relationship between 
TiO and CaH for nearby stars, Barnard's Star and KOI 
961 is illustrated in Figure|4j show ing CaH2 + CaH3 v s. 
Ti05 [see Figures 1, 3, 6 and 8 of ILepine etail (j2007)]. 
KOI 961 and Barnard's Star have similar CaH2 + CaH3 
values within their uncertainties, indicating similar effec- 
tive temperature. The position of KOI 961 to the right of 
Barnard's Star implies lower metallicity, consistent with 
the sli ghtly bluer V-K colo r. 

The ILepine et al.1 (|2007| ) TiO and CaH relations are 
diagnostics for estimating the relative differences in ef- 
fective temperature and metallicity, and useful for dis- 
tinguishing dwarfs from subdwarfs. They were not cal- 
ibrated to empirical metallicity or effective temperate 
measurements. For calibrated measurements, we used the 
if-band spectroscopic indices of RAIL 

2.3. Medium- Resolution Near-Infrared Spectroscopy 

Figure [5] shows the near- infrared spectra of KOI 961 
and Barnard's Star from the observational programs of 
Paper 1 and RAIL Both spectra were obtained with 
TripleSpec on the Palomar Observatory 200-inch Hale 
Telescope (jHerter et al.ll2008| ). We reduced the KOI 961 
spectrum with t he SpexTool program , modified for Palo- 
mar TripleSpec (|Cushing et al.l[200l ; M. Cushing, priv. 
comm.) and used spectra of nearby A0V stars for telluric 
remo val and flux calib ration with the xtellcor pack- 
age (fVacca et al.l 12003). The RAll spectrum was re- 
duced using a custom package written by P.S.M., based 
on the re duction package for the TEDI instrument, de- 
scribed in lMuirhead et al.l ()2011afj . We supplemented the 
Barnard's Star spectrum with an additional observation 
taken on UT 19 October 2011. The supplemental obser- 
vation was taken using an observing sequence allowing 
for more reliable V-band data reduction. 



TABLE 2 

Comparison of PMSU Spectroscopic Indices. 



Index — Barnard's Star — KOI 961 

1 2 

TiO 1 0.75 0.85 0.85 

TiO 2 0.64 0.63 0.75 

TiO 3 0.68 0.73 0.81 

TiO 4 0.64 0.62 0.71 

TiO 5 0.41 0.39 0.52 

CaH 1 0.74 0.84 0.77 

CaH 2 0.39 0.45 0.42 

CaH 3 0.66 0.68 0.67 

CaOH 0.37 0.45 0.54 

CTiO/CaH 3 1-37 1.87 1.26 

CxiO/CaH 4 1-43 2.02 1.33 



"Values from lReid etaD lfl995ly 

''Values calculated in this work using the spectrum from 
Dawson & Dc Robcrtis (2004). Comparing the values to 
Rcid et al. (1995), we estimate l-cr uncertainties of 0.05 for all in- 
dices. 

c The CTiO/CaH index is introduced in ILepine et alj (2007). 

''We also report the CxiO /CaH index using the coefficients from 
IDhital et aLl (120111 '). 




0.5 1 , , , i , , , i , , , i , , , i , , , i i , , 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Ti05 

Fig. 4.— PMSU spectral indices CaH2 + Ca H3 vs. Ti05 
[see Figures 1, 3, 6 and 8 of ILepine et al. (200 3)1, including all 
stars with values re ported in the PMSU survey (Rcid et aLl ll995l ; 
IHawlev et al.l [19971 ). two estimates for Barnard's Star (two black 
circles, see text), KOI 961 (red square), and t he m ean fit to the 
Solar-metallicity disk dwarfs from Lepinc et al. (2007, dashed line). 
KOI 961 has matching CaH depths but shallower Ti O depths com- 
pared to Barnard's Star. Using the interpretation ILepine et alj 
(2007), this indicates that KOI 961 has similar effective temper- 
ature but slightly lower metallicity compared to Barnard's Star, 
consistent with the i\-band measurements from RAll and Paper 
1. 

Figure [5] shows spectra for Y-, J-, H- and i-T-bands. 
The throughput calibration in Y and J bands is incom- 
plete, so we continuum normalized by fitting a quadratic 
function to the spectrum within each of the bands and 
then dividing by the fit. Regions between the bands are 
heavily obscured by telluric lines. The strong correspon- 
dence between the stars' spectral shapes and line-depths 
further reinforces the conclusion that they have similar 
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Fig. 5. — Comparison of medium-resolution near-infrared spectra 
of KOI 961 (red) and Barnard's Star (black). 1st Panel: y-band 
TripleSpec spectra of both stars. 2nd Panel: J-band TripleSpec 
spectra of both stars. The Y- and J-band spectra have each been 
continuum normalized by a 2nd degree polynomial because of in- 
complete throughput calibration. 3rd Panel: H-band spectra of 
both stars. Jith Panel: X-band spectra of both stars. The H and 
7^-band spectra of Barnard's Star were each multiplied by a nor- 
malization factor to match KOI 961. Notice the strong agreemen t 
in line depths (we identified the lines using Cushing ct al. 2005), 
and the similar spectral shape in H and K bands. The strong corre- 
spondence further reinforces the conclusion that they have similar 
effective temperatures and metallicities. 

effective temperatures and metallicities. 

2.4. High-Resolution Optical Spectroscopy 

We obtained high- resolution optical spectra of KOI 961 
with the High Resolution Echelle Spectrom eter (HIRES) 
on the Keck I Telescope (jVogt et al.lll994f) on UT 9 Oc- 
tober 2011. We exposed for 1060 seconds using settings 
that achieved a resolving power of 55000 from 3630 to 
8900 A, and a signal-to-noise ratio of 25 at 6400 A. 

Barnard's Star is currently monitored for orbiting ex- 
oplanets as part of the California Planet Search, a preci- 
sion radial- velocity survey using an iodine absorption cell 
for calibration (e.g. Uohnson et aLll2010bD . We combined 



all of the spectra taken of Barnard's Star over the course 
of the survey using regions that are not contaminated by 
molecular iodine absorption. The result is a spectrum 
with an average signal-to- noise ratio of roughly 2500. 

Figure M shows the HIRES spectra of KOI 961 and 
Barnard's Star in molecular bands corresponding to se- 
lect PMSU indices, including Ha. The correspondence 
is remarkable; however, unlike the moderate-resolution 
spectrum in Figure 02 the HIRES spectra are contin- 
uum normalized within each order using a quadratic fit. 
This will hide differences in the bulk depths of molecular 
bands, which are more apparent at moderate resolution. 
Nevertheless, the striking similarity attests to the pres- 
ence of similar molecular opacities, which are a strong 
function of abundance and temperature in M dwarfs. 
Like the mode rate-resolut i on spe ctra, neither star shows 
Ha emission. iWest et al.l (|2004l ) and iWest et all (|2008[ ) 
showed that the fraction of mid-to-late M dwarfs with 
Ha-emission decreases with height above the galactic 
plane, and the lack of Ha in either star indicates that 
they are part of an old, high disk-scale-hei ght popula- 
tion, c onsistent with their low-metallicities. We st et al.l 
(2008) estimate the activity lifetime of M4 dwarfs to be 
4.5 ±5'q Gyr, providing a lower limit on the ages of KOI 
961 and Barnard's Star. 

2.4.1. Projected Rotational Velocity 

upper limit to the projected rotational velocity of 
B arnard's star, finding V sin(i) j 2.5 km/s. 

iBrowning et al.l (|2010D report upper limit to the pro- 
jected rotational velocity of Barnard's star, finding 
Vsin(z) < 2.5 km s — 1 . This value represents the lowest 
V sin(z) measurable given the resolving power of HIRES. 
KOI 961 shows no appreciable rotational broadening 
compared to Barnard's Star, indicating similar slow rota- 
tion. To quanti fy the rotational broade ning, we followed 
the approach of [Browning et~aT1 poll . 

For each of the spectral regions in Figure [6j we 
cross-correlated the spectrum of Barnard's Star with a 
rotationally-broadened version of itself, using a convo- 
lution kernel of specified Vsin(i). We did this for sev- 
eral V sin(i) values from 0.1 to 10 km s -1 . We fitted 
Gaussian functions to the central peaks of the result- 
ing cross-correlation functions (XCFs) , and recorded the 
full-width-at-half-maximum (FWHM) of each fit. Fig- 
ure [7] shows the resulting XCF FWHM vs. V sin(i) of 
the rotational kernels for the six spectral regions (solid 
blue lines). This provides a calibration to convert XCF 
FWHM into V sin(z) for each spectral region. 

We then cross-correlated the spectrum of Barnard's 
Star with KOI 961, and interpolate that value onto the 
FWHM-ysin(z) curves for each spectral region (dashed 
red lines in Figure [7]) . The mean V sin(i) for the six re- 
gions if 2.9 km s -1 , and the standard deviation is 0.4 
km s _1 , showing that both Barnard's Star and KOI 961 
are slowly rotating. We caution that this V sin(i) value 
should be treated as an upper limit rather than an abso- 
lute measurement, since this is near the resolving power 
of HIRES, and we estimate that upper limit to be one 
sigma higher than the mean: 3.3 km s _1 . 

Assuming the inclination of the rotational axis of the 
star is 90 degress, and the 0.17 Rq stellar radius calcu- 
lated in the following section, we calculated the rotation 
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Fig, 6. — Comparison of HIRES spectra of KOI 961 (red) and Barnard's Star (black) for select molecular bands identified in Rei d et al,l 
111993 ) and H-a. The spectra are continuum normalized within each spectral region. The spectra are remarkably similar, and neither shows 
appreciable rotational broadening. The discrepancy in Ti03 (top right) is attributed to errors in the normalization of the spectra, as the 
Ti03 region is close to the edge of a HIRES order. Neithe r spectrum shows Ha emission, which has been shown to correlat e with galactic 
disk height and age for M dwarfs (West ct al. 2004, 2008); however, Barnard's Star is known for occasional Ha flaring (e.g. Pauls on et al.l 
|200H) . 



period of KOI 961 to be greater than 2.6 days. Such 
slow rotation is consistent with KOI 961 bei ng an older 
mid-M dwarf (> 1-2 Gyr, Figure 12 from llrwin et al.1 
l2011a| ): however, the low V sin(i) measurement does not 
constrain the age as well as the lack of H-a (> 4.5 Gyr), 
or the kinematic properties (see Section [3|) . 

2.4.2. Absolute Radial Velocity 

Without the iodine cell, the wavelength-solution of 
HIRES can be accurately ca librated to 0.1 km s " 1 using 
the telluric A and B bands (jChubak et al.ll2012f ). Using 
the telluric calibration, we measured the radial veloc- 
ity difference between KOI 961 and Barnard's Star by 
cross-correlating the wavelength-calibrated spectrum of 
KOI 961 with an iodine-free "template" HIRES spectrum 
of Barnard's Star taken as part of the California Planet 
Search. Applying the difference to a previously measured 
absolute radial velocity of Barnard's Star (V r = -110.85 
± 0.23 km s~ 1 . lMarcv fc Benitzlll989l ). we measured the 
absolute radial velocity of KOI 961 to be -84.48 ± 0.2 
km s . 

3. ESTIMATION OF STELLAR PARAMETERS FOR KOI 961 

The qualitative similarities between KOI 961 and 
Barnard's Star are striking. However, to accurately es- 
timate the stellar parameters of KOI 961-in particular 
mass, radius and luminosity-we must use quantitative, 
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Fig. 7. — Cross-correlation full-width-at-half-maximum (XCF 
FWHM) vs. Vsin(i). Solid blue lines: XCF FWHM for Barnard's 
Star cross-correlated with a rotationally-broadened version of itself 
for each of the six spectral regions in Figure [6] Dashed red lines: 
XCF FWHM for KOI 961 cross-correlated with Barnard's Star, 
interpolated onto the the blue lines to measure Vsin(i). 



calibrated measurements. For this purpose, we used the 
-fT-band metallicity and effective temperature measure- 
ments of the two stars from RA11 and Paper 1, which 
are consistent with the qualitative observations of Sec- 
tion [2] Table [3] summarizes the iiT-band measurements. 
In Paper 1, we interpolated the X-band measurements 
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TABLE 3 

Comparison of Measured Parameters with K-band Spectra 



Parameter 


Barnard's Star 


KOI 961 


T e ff 


3266 ± 29 K 1 


3200 ± 65 K 


[M/H] 


-0.27 ± 0.12 


-0.33 ± 0.12 


[Fe/H] 


-0.39 ± 0.17 


-0.48 ± 0.17 



"The K-hand T e g measurements a re calibrated using the BT- 
Settl-2010 model stellar atmospheres IjAllard et al.1120101 '). The un- 
certainty in T e ff is calculated using Monte Carlo simulations of 
noise in a measured spectrum. However, the uncertainty in the 
calibration itself is not currently included. 

of effective temperature and metallicity onto evolution- 
ary isochrones to estimate the stellar parameters of low- 
mass KOls, including KOI 961. However, empirical mea- 
surements of low-mass stellar parameters and predictions 
from evolutionary isochrones show significant discrepan- 
cies, particularly wit h regard to stellar radiu s and ef- 
fective temperature ([Ribasl 120061 : iTorresI 1201 If) . In this 
paper, we use evolutionary isochrones to estimate the 
differences in stellar parameters between Barnard's Star 
and KOI 961. By differencing the interpolated param- 
eters, systematic offsets in the models or measurements 
are significantly reduced. We then apply those differ- 
ences to the empirical parameters of Barnard's Star to 
estimate the parameters of KOI 961, and corresponding 
uncertainties. By "empirical," we mean measurements 
that do not depend on stellar evolutionary models. 

3.1. Summary of Empirical Measurements of Barnard's 

Star 

Unlike most known M dwarfs, Barnard's Star is bright 
enough to have a reli able parallax mea surement from 
the Hipparcos satellite (jvan Leeuw cn 2007), providing ab- 
solute photometri c magnitudes f o r the bands described 
in Section 12.11 Del fosse et al.l (|2000ft derived mass- 
luminosity relations for M dwarfs using nearby M-dwarf 
eclipsing binary systems; of their relations, the M^ CIT - 
mass relation has the lowest residuals with the least de- 
pendence on stellar metallicity (see Figure 3 from their 
paper). We calculated Mk^ by comb ining the Ken 
measurement reported in iLeggettl (|1992D with the Hip- 
parcos parallax. Using the empirical Afft- CIT -mass rela- 
tion of lDelfosse et al.l 12000), we calculated a stellar mass 
o f 0.158 M w for Barna rd's Star. 

Delf osse et al.l ([20000 do not report the systematic un- 
certainty in the M/<- CIT -mass relation. To estimate the 
uncertainty in the relation, we calculated the root-mean- 
square (RMS) of the residuals between Mk cit of the 
calibration stars and those of the relation, and found it 
to be 0.18 magnitudes. This is much larger than the un- 
certainty in Mk c1t of Barnard's Star, and will dominate 
the uncertainty in its measured mass. 

We estimate the uncertainty in the mass of Barnard's 
Star using a Monte Carlo approach. We recalculated the 
mass of Barnard's Star for 500 instances of Mk cit , each 
with noise added, drawn from a Gaussian distribution 
of standard deviation of 0.18 magnitudes. The resulting 
standard deviation of the mass distribution is 0.013 M Q . 

The radius of Barnard's Star has been measured 
previously using optical long-baseline interferometry. 
lLane et all ([2O0"lt ) reported a radius of 0.201 ± 0.008 R Q 
for Barnard's Star, using observations with the Palomar 



Testbed Interferometer. The uncertainty in the stellar 
radius comes from a combination of measurement errors, 
uncertainty in the angular size of the calibrator stars and 
uncertainties in t he ado pted limb- darkening parameters. 
iSegransan et al.l ([20031 ) recalculated the radius using a 
different set of limb-darkening coefficients, and reported 
a radius of 0.196 ± . 008 R q for Barnard's Star, consis- 
tent with lLane et al.l ((2001) to within the quoted uncer- 
tainty. This shows that the interferometrically measured 
stellar diameter is robust against assumptions about the 
limb darkening and strengthens our confidence that this 
radius measurement represents an accurate measurement 

of the tru e radius of Barnard's sta r. 

Finally, iDawson &: De Robertisl (|2004f ) combined and 
calibrated available spectra of Barnard's Star across op- 
tical and infrared wavelengths to accurately measure its 
luminosity [(3.46 ± 0.17) x 10 -3 L Q ], which, when com- 
bined with the radius measurement, provides an empiri- 
cal measurement of the star's effective temperature (3134 
± 102 K). This agrees remarkably well with the effective 
temperature measured using X-band in RA11 of 3266 ± 
29 K. The quoted uncertainty for the effective tempera- 
ture measured in RA11 only accounts for measurement 
errors in the i^-band spectrum, and not systematic ef- 
fects from errors in the calibration. 

3.2. Perturbation Analysis 

To estimate the difference in mass, radius and lumi- 
nosity between Barnard's Star and KOI 961, we inter- 
polated the .fT-band measurements of stellar effective 
te mperature and metall icity onto the 5-G yr isochrones 
of IBaraffe et al.l (|1998f ). We use the (Baraffe et all 
(1998) models because, unlike other sets of evolution- 
ary isochrones, they include grid points with stellar ef- 
fective tem peratures below 320 K for sub-solar metal- 
licity. The IBaraffe et al.l (| 19981 ) isochrones are available 
in two metallicities, [M/H] = 0.0 and -0.5, which bracket 
the measured metallicities of Barnard's Star and KOI 
961. The [M/H] = -0.5 isochrones are only available for 
a mixing-length parameter (a) of 1 and for a Helium 
abundance ([Y]) of 0.25, so we choose the [M/H] = 0.0 
isochrones with those same parameters. To test the ef- 
fect of age on the parameter estimation, we repeated the 
calculations using 10-Gyr isochrones. The resulting stel- 
lar parameters were only marginally different compared 
to the uncertainties, a consequence of the slow evolution 
of M dwarfs. 

Figure [8] shows stellar radius and stellar mass 
vs. effective tempe rature for the 5-Gyr isochrones of 
IBaraffe et al.l (|1998f) and the interpolated location of 
Barnard's Star and KOI 961 based on the if-band mea- 
surements. We applied the difference in mass, radius and 
luminosity between the model interpolants to the em- 
pirical measurements of Barnard's Star to estimate the 
values for KOI 961. 

To estimate the uncertainties in KOI 961 's parameters, 
we interpolated onto the isochrones 500 Monte Carlo sim- 
ulations of the if-band T c ff and [M/H] measurements, 
each with random noise added. The noise is drawn from 
a normal distribution with standard deviation of 100 K 
for T cf f and 0.12 dex for [M/H]. We then calculated the 
standard deviation of the distribution in the resultant 
radii, masses and luminosities. We chose to inflate the 
measurement uncertainties in T c g to 100 K to account 
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Fig. 8. — Top Panel: Stellar radius vs. effective temperature. 
Bottom Panel: Stellar mass vs. effective temperature. Colors de- 
note metallity , [M/H l. Lines: The 5-Gyr evolutionary isochrones of 
IBaraffe et aU(Tl998T) for [M/H] = 0.0 and [M/H] = -0.5. Large col- 
ored circle and square: Interpolated position of Barnard's Stars and 
KOI 961, respectively, based on if -band measurements of metal- 
licity and effective temperature. Small colored circles and squares: 
Monte Carlo simulations of noise in the i\-band measurements, 
used to estimate the error in the mass, radius and luminosity of the 
interpolated values. Black circle: Empirically measured values for 
Barnard's Star and uncertainties. Black square: Our estimate for 
the radius and effective temperature of KOI 961, found by applying 
the difference in interpolated values to the empirical measurements 
of Barnard's Star. 



for potential errors in the calibration. For [M/H] we 
used an uncertainty of 0.12 based on the calibrat ion un- 
certainty calculated in IRoias-Avala et ail ()2010D . This 
provides uncertainties in the interpolated masses, radii 
and luminosities of KOI 961 and Barnard's Star, which, 
when added in quadrature, provide the uncertainty in 
the difference between the stars' parameters. We then 
added in quadrature the uncertainties in the differences 
to the empirical measurements of Barnard's Star, which 
provides the uncertainty of KOI 961's parameters. Table 
0] lists the empirical measurements of Barnard's Star and 
our estimated values for the stellar parameters of KOI 
961 and their uncertainties. 

We calculated the space motion of KO I 961 by combin- 
ing th e proper motion measurements of ILepine fc Shara 
(2005) with our estimated distance and measured radial 
velocity. Figure [§] shows a Toomre diagram of nearby 
G- K- and M-type primary stars with measured trigono- 
metric parallaxes gre ater than 100 mas . Using the pop- 
ulation definitions of iFuhrmannl (|2004D . it is clear that 
both Barnard's Star and KOI 961 have kinematic prop- 
erties consistent with the "thick disk" , containing stars 
with ages of 7 to 13 Gyr dFeltzing fc Bensbvi 12008). We 
acknowledge that the presence of a thick disk in the 
Milky Way is controversial, however the difference in 



TABLE 4 

Barnard's Star and KOI 961 Stellar Parameters 



Parameter 


Barnard's Star 
(empirically measured) 


KOI 961 
(this work) 


R* 


0.199 ± (J.UUo Hq 


(J. 17 ± U.U4 Hq 


M 


n 1 ^ -i- n m ^ 




L* 


(3.46 ± 0.17) x 10" 3 L Q 3 


2.40 x l0- 3 0±0 - 3 L 4 


T e ff 


3134 ± 102 K 3 


3068 ± 174 K 


Mk, cit 


8.21 ± 0.03 5 ' 6 


8.55 ± 0.35 


M Ks 


8.22 ± 0.03 5 - 7 


8.53 ± 0.35 


d 


1.824 ± 0.005 pc 7 


38.7 ± 6.3 pc 


U 


-134 i 1km s" 1 


50 ±10 km s" 1 


V 


18 ± 1 km s" 1 


-70 ± 1 km s" 1 


w 


27 ± 1 km s" 1 


-59 ± 8 km s" 1 


V sin(j) 


< 2.5 km s~ 18 


< 3.3 km s" 1 



"Weighed mean of lLane et alj (200li 0.201 ± 0.008 Rq) and 
ISeeransan eTaLl J200J 0.196 ± 0.008 Rq) ^ 

'Calculated using the Mx,ciT _mass relations of IDelfosse et alj 
1200TJ). 

IDawson fc De Robert^ 1 20041) 

d Since luminosity is a strong power of the effective temperature 
and radius, the distribution of interpolated Monte Carlo luminosi- 
ties is asymmetric. Therefore, we quote the error in the inferred 
luminosity of KOI 961 in the exponent, equivalent to the standard 
deviation of the log of the res u ltant luminosity distribution. 

^Parallax from van Lccuwcn ( 2 0071) 

jLeggetti g99g) 

»2MASS, Cutri eTaLl J2003D 

' Browning et al. (2010) 
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Fig. 9. — Toomre Diagram of G- K- and M-type primary stars 
with measured trigonometric parallaxes > 100 mas. We i nclud e 
the thin and thick disk boundaries defined by Fuhrmann 1120041 ) . 
Barnard's Star (red circle) and KOI 961 (blue square) are indicated 
and have space velocities consistent with thic k-disk stars. We cor- 
rected for the solar motion using the values of Francis & Anderson 
(2009): U = 7.5 + -1.0, V = 13.5 + 0.3, W = 6.8 + -0.1 km/s. 



stellar kinematics as a function of di sk scale-height th e 
galactic disk is firmly established (e.g. lBovv et al.ll2011|) . 
The kinematic properties, metallicities, Ha emission and 
Vsin(i) measurements of KOI 961 and Barnard's Star 
are all consistent with an old, high disk-scale-height stel- 
lar population. 



3.3. Consistency Check on GJ 1214 

As a consistency check of our methods, we performed 
an identical perturbation analysis on GJ 1214, a low- 
mass star with a transiting Super-Earth: GJ 1214 b 
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(jCharbonneau et al.l [2009h . iCarter et all (j201ibD mea- 
sured multiple transit light-curves of GJ 1214 b and cal- 
culated the mean de nsity of the star using the fitted tran- 
sit parameters (e.g. iSeager & : Mallcn-Ornclas 200 3|). GJ 
1214 has a parallax measurement by Ivan Altena et aLl 
(|1995[) and 2MASS J HK photometry Com bining the 
mass-Mx relation of iDelfosse et al.l (120001) with the 
transit-derived stellar density. ICarter et al.l ()2011bf ) re- 
ported a stellar mass of 0.157 ± 0.012 M Q and a stel- 
lar radius of 0.210 ± 0.007 R Q for GJ 1214 (Method 
A from their paper). The stellar radius was larger 
than estimated usin g photometric relationships from the 
iBaraffe eT al. (1998) evolutionary isochrones, which they 
also report as i?* = 0.179 ± 0.006 R© (Method B from 
their paper). 

RAll report a stellar effective temperature of 3245 
± 31 K and metallicity ([M/H]) of 0.15 ± 0.12 for GJ 
1214 using the if -band spectroscopic indices. With an 
effective temperature similar to Barnard's Star and KOI 
961, as well as available if-band measurements, GJ 1214 
provides a valuable opportunity to test our perturbation 
analysis. 15 

We proceed assuming Solar-metallicity for GJ 1214, 
since we cannot interpolate the super-Solar if-band value 
of_[M/H] = +0.15 onto the isochrones of IBaraffe et al.l 
(1998). Using the same perturbation procedure used for 
KOI 961, and the if -band effective temperature mea- 
surement, we calculate a a stellar mass of 0.17 ± 0.05 
Mq and a stellar radius of 0.21 ± 0.04 R for GJ 1214. 
The values are remarkably c o nsisten t with the results 
of Method A by ICarter et all (|2011bl ) to within our es- 
timate uncertainties. Incorporating the if -band metal- 
licity measurement of GJ 1214 would result in a larger 
radius had we used super-Solar metallicity isochrones, 
but not by greater than Ai?* = 0.1 R (see Figure [8]): 
still well within the estimated uncertainty. We conclude 
that the perturbation analysis is an effective method for 
determining accurate stellar parameters of M dwarfs in 
the absence of parallax measurements. 

4. TRANSIT PARAMETERS 

Now we turn to a discussion o f the t ransits detected 
around KOI 961. iBorucki et al.l (|2011l) reported tran- 
sit parameters (period, duration, impact parameter and 
planet-star radius ratio) for the three planet-candidates 
orbiting KOI 961. However, the published radius and 
impact parameter for KOI 961.02 are inconsistent: the 
impact parameter and planet-to-star radius ratio in- 
dicate that the transit should not occur (6/i?*=1.29, 
i?p/i?*=0.194). With this in mind, we chose to per- 
form an independent fit to the public Kepler data for the 
three planet-candidates. In the following section, we use 

15 Other possible test stars include the short-period M dwarf 
eclipsing-binaries CM Draconis and KOI 126, which have sim- 
ilar effective temperatures to Barnard's Star and KOI 961 (~ 
3000 K) and highly accurate mass and radius measurements 
IILacv|ll977l : IMetcalfe eT"aLlll999 : IMorales et~aTll2009l ; ICarter et all 
l2011aT ). However, it has been suggested that mutual magnetic in- 
teractions in short-period M dwarf binaries introduce discrepancies 
in the mass-radius-effecti ve temperature relatio nships as compared 
to single field M dwarfs (Chabricr et al. 2 0071') ■ and recent obser- 
vational evidence supports this scenario (Kraus et al. 2011). Also, 
measuring the A"- band T c g and [M/H] indices of individual stars in 
unresolved, short-period eclipsing binaries is challenging and has 
not yet been done. Therefore, we believe GJ 1214 is the best test 
object. 



the resulting parameters to calculate the probability that 
the transit signals are not planetary in nature, and we 
find this probability is less than 1% for all three signals. 
Therefore, in this section, we refer to them as planets 
rather than candidates. 

We use Levenberg-Marquardt least-square minimiza- 
tion to refit the light curve of KOI 961 and determine re- 
vised transit parameters for all three planets. This task is 
complicated by a number of simultaneous transits of two 
or more planets. For these situations, we assume that the 
total loss of light is the superposition of the losses com- 
puted individually for each planet. We do not attempt 
to model the short brightening in the light curve as a 
result of the possible m utual alignment of three or m ore 
bodies [as described by (|Ragozzine fc Holmad[2010l )]. 

In principle, one could describe the transit model for 
each planet with three minimal parameters: the planet 
to star radius ratio, the orbital inclination, and the in- 
stantaneous separation of planet and star, normalized to 
the stellar radius: a/R*. Additionally, we may specify 
the ephemeris and period for each planet and parame- 
ters describing the radial brightness profile of the star. 
In the case of KOI-961, each individual transit duration 
is comparable to the integration time for the long cadence 
Kepler data. As a result, there exists a near exact de- 
generacy between these three minimal parameters. For- 
tunately, the normalized instantaneous separation can be 
constrained given our knowledge of the star if we are will- 
ing to assume the orbits of all three planets are circular. 
In this case, a/i?* is specified in terms of the mean stel- 
lar density and the orbital period. For our fitted transit 
models, we assume circular orbits and apply the mean 
stellar density calculated from the stellar mass and ra- 
dius determinations in Section [3] 

We report in Table [4] the best-fitting radius ratios and 
their errors along with the best-fitting impact parame- 
ter subject to the assumption of circular orbits and the 
specified mean stellar density. In Figure [TU1 we show the 
Kepler light curve for each planet after folding the data 
on the best-fitting linear ephemeris and after having re- 
moved the transit profiles of the remaining two planets. 
We also plot representative best-fitting transit models - 
both the convolved and the continuous profile - for a few 
impact parameters. The corresponding, best-fitting ra- 
dius ratios for each impact parameter are listed in each 
panel. High im pact parameter solut ions (> 1), as listed 
for KOI-961 in IBorucki et al.l ()2011l ). are disfavored by 
the data under these constraints for all three planets. As 
a result, we determine a significantly smaller radius ratio 
for all three planets than previously reported. 

Combining the best fit transit parameters with our new 
stellar parameters, we report the planets' physical radii, 
equilibrium temperatures and corresponding uncertain- 
ties in Table |U All three planets have radii smaller than 
the Earth's, and KOI 961.03 is 1.07 ± 0.33 times the 
radius of Mars. The uncertainties in the planet radii 
are dominated by the uncertainties in the stellar radius. 
To calculate the equilibrium temperatures of the planets, 
we scaled the equilibrium temperature of Earth at 255 
K-which assumes a 30% albedo and 100% re-radiation 
fraction-to the semi-major axes of the planetary orbits 
and the luminosity of KOI 961 and propagated the un- 
certainties in those parameters. 
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TABLE 5 

KOI 961 Transit Parameters and Planet Parameters 



KOI 


t [BJD - 2454900] 


Period [days] 


a [AU] 


Rp/R* 


b [R*] 


Rp [Re] 


T cq [K] 


961.01 


103.48329 ± 5.4xl0~ 4 


1.2137672 ± 4.6 xlO" 6 


0.0116 


0.0419 ± 0.0018 


0.711 ± 0.049 


0.78 ± 0.22 


519 ± 52 


961.02 


66.86952 ± 4.2xl0~ 4 


0.45328509 ± 9.7 xl0~ 7 


0.0060 


0.0395 ± 0.0011 


0.520 ± 0.068 


0.73 ± 0.20 


720 ± 73 


961.03 


66.7847 ± 1.3xl0" 3 


1.865169 ± 1.4xl0" 5 


0.0154 


0.0307 ± 0.0019 


0.681 ± 0.085 


0.57 ± 0.18 


450 ± 45 



1.001 
1.000 
0.999 
0.998 
0.997 

1.001 
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Fig. 10. — Kepler light curves for each planet after folding the 
data on the best-fitting linear ephemeris and after having removed 
the best-fitting transit profiles of the remaining two planets. We in- 
clude representative best-fitting light curves - both the convolved 
and the continuous profile (solid lines and dotted lines, respec- 
tively) - for three impact parameters: (red), 0.6 (green) and 0.9 
(blue). The corresponding, best-fitting radius ratios for each im- 
pact parameter are listed in each panel . High impact param eter 
solutions (> 1), as listed for KOI-961 in Borucki ct al. (2011), are 
disfavored by the data under these constraints, evidenced by the 
poor match between the solid blue curve and the Kepler data at 
the start and end of each folded light curve. 



5. THE PLANETARY NATURE OF THE TRANSITS 

It is conventional to call transit signals candidates until 
the planets are dynamically confirmed, either by radial 
velocity or transit timing measurements. However, we 
do not have Doppler spectroscopy for KOI 961 and the 
transit signals show no obvious timing variations. For- 
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Fig. 11. — The measured fitting statistic, x 2 i versus the mod- 
eled impact parameter for each planet (black). In all cases large 
impact parameters are disfavored by the iCepler photometry, and 
in the case of planets .01 and .03, the light curves provide strong 
lower constraints on the impact parameter. We include the best- 
fitting planet-to-star radius ratio for each impact parameter (red), 
illustrating the small spread in radius ratios across the minimum 
of the x 2 vs - impact parameter curve. The best-fitting impact 
parameters and radius ratios are listed in Table [4] 



tunatcly, given that there are three transit signals, each 
with a period of less than 2 days, it would require an 
extremely contrived scenario to account for all of them 
by non-planetary means. In the following section we 
quantify the probability that each of the signals may be 
an astrophysical false positive, demonstrating that this 
probability is low enough to consider all three planets 
validated. We also discuss the constraints on potential 
diluting third light in the aperture, which would lead to 
underestimating the planetary radii. 

5.1. Observational constraints on false positive 
scenarios 

Photometric signals that appear to be transiting plan- 
ets may be caused by astrophysical scenarios other than 
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a planet orbiting the target star. In this work, we con- 
sider a false positive to be any scenario different from 
a planet transiting a star in the KOI 961 system. This 
may be either an eclipsing binary or a transiting planet 
around a chance-aligned background star or an eclipsing 
binary that is physically bound to KOI 961 in a hierar- 
chical triple system. Each of these scenarios involves the 
presence of another star within the Kepler aperture of 
KOI 961. 

We constrain the properties of such a star using 
archival and recently acquired imaging data, which we 
show in Figure [T2] The large proper motion of KOI 961 
allows us to use digitized ima ges from the Palomar Ob- 
servatory Sky S urvey (POSS; (Minkowski fe AbeHlll963t 
Rci d et al.lll991l ) to rule out the presence of any contami- 
nating object up to 4.8 magnitudes fainter than KOI 961 
in the blue-sensitive plates, anywhere with in the Kepler 
apertures used in the iBorucki et al.1 ([20111 ) data release. 
We also use both publ icly available J-band UKIRT data 
(|Lawrence et al.ll2007|) and a HIRES guider camera im- 
age we obtained using Keck I (equivalent of combined R- 
and U-band response), which allows us to constrain the 
presence of blends within a few arcseconds of the current 
position of KOI 961. 

In addition, the HIRES spectrum of KOI 961 con- 
strains the presence of any potentially blending com- 
panion through the absence of noticeably contaminating 
lines. The symmetry of the cross-correlation profile of the 
spectrum with that of Barnard's Star reveals that KOI 
961 has no companion of comparable brightness (AV~1; 
or mass ratio q ~ 0.75) with velocity separation >2 km 
s" 1 re.g. lHoward"et^l2T)l7ih. 16 

Using the methods of lMorton fe Johnson! (j 2 1 lh , which 
use simulations of Galactic structure and stellar evolu- 
tion to estimate the sky density of background stars, 
we determine that there is a probability of about 0.13 
for there to be a chance- aligned star undetected by the 
POSS images or the Keck guider i mage anywhere within 
the Kepler apertures used for the IBorucki et al.l ([20111 ) 
data release. To determine the extent to which the above 
observations constrain hierarchical triple scenarios, we 
simulate a population of stellar binary companions with 
mass ratio, period, and eccentri city distributions typi- 
cal of local multiplici ty surveys (jRaghavan et al.1 [20101 : 
iTokovinin et al.l I2010T ) . and assign magnitudes t o eac h 
companion based on the models of iBaraffe et al.l (|1998j ). 
We then place each simulated companion at a random 
point on its orbit (using a uniform distribution of mean 
anomaly) and use the resulting projected separation (as- 
suming a distance of 40 pc) , relative velocity, and mag- 
nitude difference to determine whether it would be ob- 
servable given our observations. We find that 59% of 
companions to KOI 961 would have eluded detection. As- 
suming an overall 40% binary fraction for M dwarfs (e.g. 
iFischer fe Marcvlll992t iCovev et all 120081 ). this yields a 
probability of 0.24 that KOI 961 has a undetected com- 
panion. 

5.2. False Positive Probability 

16 Incidentally, this observational constraint also eliminates the 
false positive scenario (not discussed in the text) in which one of 
the transit signals is caused by a similar-mass binary companion 
to KOI 961 with a grazing eclipse. 




(a) POSS-I blue, 1951 



(b) POSS-I red, 1951 
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(d) POSS-II red, 1991 




(e) HIRES Slit Viewer (this (f) HIRES Slit Viewer close- 
work), 2011 up of KOI 961 

Fig. 12. — Images of KOI 961 from the Palomar Observatory Sky 
Survey (POSS, p anels a, b, c and d, e.g. [Minkowski & AbeTH 119631 ; 
IReid e t al. 199lf) and from the HIRES slit -viewer (this work, pan- 
els e and f). Panels a through e are roughly 1 arcminute on each 
side and show the large proper motion of KOI 961 over 60 years. 
The magenta circles indicate the current (2011) position of KOI 
961 and have diam eters of 6", roughly equal to the poin t-spread- 
function of Kepler l IBrvson et al.ll201ot [Koch et al.ll2010D . We also 
include the apertures used by Kepler for flux extraction of KOI 
961 during quarters 1, 2 and 3, which we obtained from the Mul- 
timission Archive at the Space Telescope Science Institute [green 
outlines). We rule out the presence of a contaminating star < 4.8 
magnitudes fainter than KOI 961 using the POSS 1 blue-sensitive 
image. We use the POSS images and HIRES slit-viewer image to 
limit potential false-positive scenarios for the transit signals. 



Calculating the false positive probability (FPP) for a 
particular transit signal requires comparing the proba- 
bilities of astrophysical false positive scenarios with an 
assumed probability of the planet hypothesis. To this 
end, we calculate the probability that each of the po- 
tentially blending scenarios discussed above might be an 
eclipsing system mimicking a planet-like transit signal 
around KOI 961. 

Direct applicat ion of the exact results of 

iMorton fc Johnsonl (|201lD is not appropriate for 
these signals, because the transit shape does not imme- 
diately rule out grazing eclipses. As we show in Section 
ITOl the 29.4-minute observing cadence of the public 
Kepler data strongly affects the apparent shape of the 
light curve, including its dep th (see Figure ITUl). I n stead , 
we extend the methods of IMorton fc Johnsonl ([201 ID 
tailored to the specifics of the KOI 961 system. 
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TABLE 6 

False Positive Probabilities for KOI 961 





BGEB 1 


HEB 2 


BGpl 3 


FPP 4 


961.01 
961.02 
961.03 


9.7e-05 
2.8e-04 
l.le-03 


3.8e-05 
1.5e-05 
1.2e-03 


5.4e-04 
l.lc-03 
l.le-03 


6.7c-04 
1.4e-03 
3.4e-03 



"Background eclipsing binary 
'Hierarchical eclipsing binary 
c Transiting planet around a background star 
d Total false positive probability 
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Fig. 13. — The probability density as a function of observed transit duration and depth for realistic populations of the three astrophysical 
false positive scenarios discussed in the text a s well as the bona fide planet scenario. We choose a statistically representative sample of 
configurations and model the light curves using Mandcl & Agol (2002), convolved with a 29.4 minute cadence, assuming all circular orbits. 
Included in these density plots are only those scenarios in which secondary eclipses would not be detectable in the KOI 961 light curve. 
The observed transit signals of KOIs 961.01, 961.02, and 961.03 are marked (with durations scaled to match the period of 961.01, which 
was the basis for the population simulations), and do not fall within a high-likelihood region for any of the false positive scenarios. The 
relative probability for each planet to be caused by each scenario is the product of the prior and likelihood for that scenario. 

uated at the specific depth and duration of the observed 
signal (the "likelihood" ) , and 7ij is the a priori probabil- 
ity of each scenario existing. 

The a priori probability for each scenario has two com- 
ponents: first, that the particular blending star config- 
uration exists; second, that it eclipses and can mimic 
a planet around KOI 961. For the background eclips- 
ing binary (BGEB) scenario, we assume that a chance- 
aligned star has a 50% chance of being a binary, and a 
10% chance of that binary being on a short-period or- 



The FPP for any of the three signals may then be writ- 
ten as follows: 

FPP = 1 _ *Pi4>i(*,r) 



where the index i represents each of the false positive sce- 
narios and the planet hypothesis, Ct(5,T) is the proba- 
bility density function (PDF) for the given scenario eval- 
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bit ([Raghavan et al.lfeOlOft . using < 50 days as a defini- 
tion of "short-period" to allow convenient comparison to 
planet occurrence studies. Using 0.18 as a conservative 
average eclipse probability (corresponding to the graz- 
ing eclipse probability of two 0.5 Rq stars orbiting each 
other with a 1.2 day period), this gives a probability of 
0.13 x 0.50 x 0.1 x 0.18 = 1.1 x 10" 3 that there is a BGEB 
within the Kepler apertures of KOI 961. In addition, the 
binary must not show a secondary eclipse in the diluted 
light curve. We determine by inspection that a signal 
with a depth of 5 x 10~ 4 would have been detectable by 
Kepler, and so we use that as a conservative limit for 
the maximum depth of any putative diluted secondary 
eclipse. We find that this disqualifies 39% of BGEBs, 
resulting in 7Tbgeb = 6.9 x 10~ 4 . 

To estimate the prior probability of a hierarchical triple 
eclipsing binary (HEB) system, we assume that a poten- 
tial companion to KOI 961 has a 40% chance of hav- 
ing a companion of its own, a 20% chance of that com- 
panion being on a short orbit [we conservatively use 
a higher short-orbit probability for a hierarchical sys- 
tem than for an isolated binary motivated by the ob- 
servation thj^jnairy_close_binaries have distant com- 
panions ([Tokovinin et al .11200 61)], and an eclipse proba- 
bility of 0.10 (corresponding to two 0.15 R© stars in a 
1.2 day orbit). Furthermore, we simulate a population 
of HEBs according to a flat mass ratio distribution and 
find that only 2.3% of HEB scenarios have diluted sec- 
ondary eclipses shallower than 5 x 10~ 4 . This results in 
ttheb = 0.24 x 0.40 x 0.20 x 0.10 x 0.023 = 4.4 x 10~ 5 . 

For the background transiting planet hypothesis 
(BGpl), we assume an overall 40% occurrence rate for 
planets with periods < 50 days ([Howard et al.l 120101 
[20T1 IBonfils et alJl20ll . and an average transit prob- 
ability of 0.13 (corresponding to a planet orbiting a 0.5 
Rq star in a 1.2 day period), giving 7rBGpi = 6.8 x 10~ 3 . 
These same assumptions (except for a transit probability 
of 0.07 corresponding to the radius of KOI 961) yield an 
a priori probability for the bona fide KOI 961 transiting 
planet hypothesis of vr p i = 0.40 x 0.07 = 0.028. We note 
that we have assumed isotropically oriented orbits for all 
the eclipse and transit probabilities we have used. This 
may not necessarily be the case for multiple planet sys- 
tems (orbits may be coplanar) or for hierarchical systems 
(orbits of a hierarchical binary may be aligned with the 
plane of the planets) , but assuming isotropic orientation 
for all orbital planes is a conservative approach. 

To calculate the likelihood factors for each of the sce- 
narios, we sim ulate the expect e d po pulation of light 
curves following iMandel fe Agoll (|2002f ) , allowing for all 
impact parameters (including grazing eclipses) , including 
limb darkening, and accounting for the effects of the 29.4- 
minute obs erving cadenc e . In contrast to the B LENDER 
procedure ([Torres! l20TTt iFressin et all 120 11 al ibi), we do 
not fit false positive light curves to the data; instead 
we choose the scenarios to simulate according to their a 
priori probability of ex istence, which we determ ine us- 
ing the methodology of lMorton fc Johnson! ([201 ll ). This 
allows us to directly extract the joint PDF of duration 
and depth for each scenario, taking care to include only 
those light curves with undetectable diluted secondary 
eclipses. We also create the same distribution for the 
planet scenarios, assuming a population of planets with 



a radius distribution dN/dR p ~ R 2 from 0.4 to 20 
and accounting for the mass and radius errors of 
KOI 961. The resu l t is ex actly analogous to Figure 5 of 
iMorton fc Johnson! ([201 ll ) extended into the duration di- 
mension. "Duration" and "depth" for our purposes here 
are determined by inspection, not by any fitting proce- 
dure; "duration" is the total time the transit model flux 
spends below unity, and "depth" is simply the minimum 
value of the model light curve during primary transit. 
These values correspond to what one would naively call 
duration and depth from inspecting a real light curve by 
eye; they have no direct physical significance except for 
providing a convenient, model-independent way to pa- 
rameterize a light curve. We then evaluate the PDFs at 
the duration and depth of each of the three candidate 
signals (integrating over errors in both dimensions) to 
obtain each of the likelihood factors (i.e. the "probabil- 
ity of the data given the model"). The simulations we 
use to create the PDFs are all based on the period of 
KOI 961.01, so in order to evaluate the likelihoods for 
961.02 and 961.03 at the proper locations we scale their 
true durations by a factor of (P/l^ld) 1 ' 3 . Figure [TBI 
illustrates each of these PDFs and shows the likelihood 
value for each candidate for each scenario. 

Table [6] lists the normalized probabilities 
(■KiCi I J2 TTiCi) for each of the false positive sce- 
narios for each of the candidates as well as the overall 
FPP. The total false positive probabilities for 961.01, 
961.02, and 961.03 are 6.7 x 10~ 4 (~1 in 1500), 
1.4 x 10~ 3 (-1 in 700), and 3.4 x 10~ 3 (-1 in 290), 
respectively. We note that these calculations have used 
conservative assumptions at every stage. In particular 
we have ignored the fact that the ensemble evidence 
of the Kepler planet sample points to coplanarity 
being co mmon among close-in , tightly packed planetary 
systems ([Lissauer et al.|[20Tll) ; if we properly accounted 
for this effect, the planet prior for each of the candidates 
would increase (corresponding to an increased transit 
probability), leading to correspondingly lower FPPs. We 
therefore consider all three candidates in the KOI 961 
system to be statistically validated planets, according to 
the standards set in the literature by the Kepler team. 

5.3. Misclassified planet probability 

We note that even though we have demonstrated all 
three transit signals in the KOI 961 system to be of 
planetary origin, there is still some probability of their 
radii being underestimated due to the presence of dilut- 
ing flux within the aperture. Specifically, there are three 
sub-scenarios that are encompassed by the "planet in the 
KOI 961 system" scenario: 

(a) KOI 961. Ox orbits KOI 961, which is a single star. 

(b) KOI 961. Ox orbits the primary component of KOI 
961, which is an unresolved binary star system. 

(c) KOI 961. Ox orbits the secondary component of KOI 
961, which is an unresolved binary star system. 

We calculated in the previous subsection that there was 
a 24% chance of KOI 961 having an undetected binary 
companion, giving scenario (a) a probability of 0.76. We 
then assume that if KOI 961 is a binary, there is equal 
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probability for a planet to be orbiting either component, 
giving both (b) and (c) probabilities of 0.12. Using the 
same binary companion simulations discussed above, we 
calculate the effect this has on the radii of the planets, 
according to the following: 



-^p,true — ^P, inferred (l 4~ 10 P ) 



(2) 



where i? D 

.true is the true planet radius, -Rp^infcrrcd is the 
radius inferred assuming the star is single, and AMkc P 
is the magnitude difference between the two components 
(which would be negative if the planet is orbiting the 
fainter star). As each binary companion in our simula- 
tion has an assigned Kepler magnitude, we can turn this 
putative population into a probability distribution func- 
tion for the true planet radius, which we plot in Figure 
[14] for KOI 961.01. 73% of the probability mass remains 
below 1 i?0 , compared to 84% of probability mass for the 
(a) scenario probability distribution alone. The planets 
961.02 and 961.03 have analogous probability distribu- 
tions, with lower mean values, so we conclude that the 
possibility of a binary companion does not alter our con- 
clusions of the sub-Earth-sized nature of the KOI 961 
planets. The resulting uncertainties in the planet radii 
due to mischaracterization are similar to the quoted un- 
certainties in Table 0] which were calculated assuming 
scenario (a) and dominated by errors in the stellar pa- 
rameters. Therefore, we proceed assuming the uncertain- 
ties in Table H 

6. MASSES, DYNAMICS AND ARCHITECTURE 

We do not currently have measurements of the plan- 
ets' masses, as KOI 961 is faint (V — 16.1 ± 0.055) and 
the expected radial-velocity semi- amplitudes are small: 
K = 1.1, 1.25 and 0.33 m s" 1 for KOI 961.01, 961.02 
and 961.03 respectively. In Figure [T5] we show a plan- 
etary mass-radius diagram with theoretical predictions 
for cold spheres composed of pure iron, r ock (Mg2Si04) 
and w ater-ice, interpolated from Table 1 of Fo rtnev et al.1 
(|2007h . and we show potential masses for KOI 961 's three 
planets. Planetary bodies arc unlikely to be dominated 
by elements that are heavier than iron; therefore, the 
pure iron mass-radius relation provides firm upper limits 
to the planet masses. We report la upper limits of 2.73, 
2.06 and 0.90 M ffi for KOI 961.01, .02 and .03 respec- 
tively. With such low-masses, we believe it is unlikely 
that the planets could retain gaseous atmospheres with 
scale-heights that are a significant fraction of the planets' 
opaque radii. The range of possible gaseous scale-heights 
for a specified radius, surface gravity and equilibrium 
temperature provide a lower-limit on the planets masses; 
however, we leave this calculation for a future paper. 

To consider the dynamical effects in this multiple- 
planet system, we assume a stellar mass 0.14 M , and 
the 1-er upper limits on planetary masses. We considered 
coplanar orbits, as this is highly likely due to the multi- 
transiting nature of this system (Lissauer et al. 2011). 
In the following calculations we also assume near circular 
orbits. The eccentricity damping times cale due to tides 
raised on the planets due to the star is (|Wu &; Goldreichl 
2002]) : 



(3) 



where Q' is a parameter th at is < 10 3 for terrestrial 
bodies (jMardling & Lin 2004|). Using measured values 
for P, a, AT*, and i?*, and assuming a terrestrial density 
to estimate M p , all the planets have short eccentricity- 
damping timescales compared to the age of the system (> 

4.5 Gyr based on H-a), even the outermost planet (r e ~ 
0.1 Gyr). This justifies the assumption of circular orbits 
in modeling the present-day dynamics of the system. 

Given the very low-mass planets we have inferred, and 
the distances from strong mean-motion resonances, we 
would not expect large transit timing signatures. We per- 
formed 3-body numerical integrations t o determine simu - 
lated transit times using the method of lFabrvckvl (|2010ft . 
Each simulation had two planets - their perturbations 
were considered pairwise. The strongest interaction was 
between KOI-961.01 and KOI-961.03, with peak-to-peak 
deviations from a linear ephemeris of 0.7 minutes and 

2.6 minutes, respectively. Other pairs showed less than 
1 second of deviation. Although KOI-961.03's deviation 
may be detected in the full Kepler data set, we expect 
other interactions to be undetectable. 

We also investigated the long-term stability of this sys- 
tem. First, we compute the separ ation in number of 
Hill spheres ((Chamber s et al.l I1996D between the inner 
two planets and between the outer two plan ets: they 
are 19 .5 and 9.6, respectively. According to iGladmanl 
pairs of planets are Hill stable (i.e. orbits that 
are initially circular may not cross) if this value is greater 
than 3.46. When three-planet systems are considered, 
the stabi lity boundary moves t o larger separation; for 
instance. iChambers et al.l (|1996t )[fig. 3] shows instability 
is possible in ~ 10 6 orbits if the Hill separation between 
pairs of planets in three-planet systems is 6. For our 
case, scaling relations suggest that stabi lity should easily 
be ac hieved in th is three-planet system (|Chambers et al.l 
1996). Moreover. iLissauer et al.l ([201 1[ ) integrated three- 
planet systems with precession mimicking the general 
relativistic effects, including a model for KOI-961 with 
much larger mass ratios than we suggest. They found 
stability for over 10 billion orbits. Future work could 
consider precession due to tidal and rotational effects and 
initially non-circular orbits. 

7. SUMMARY AND DISCUSSION 

The Kepler Mission has demonstrated unprecedented 
photometric precision over long time baselines, opening 
up opportunities to study planets with r adii compara- 
ble to the Solar System te rrestrial planets (jBorucki et al.l 
l2011tlHoward et al.ll201lD . However, a transit light curve 
does not provide a direct measure of a planets radius, but 
instead gives only the ratio of planet and stellar radii. 
Thus, the precision of the Kepler transit light curves 
cannot be fully exploited without precise and accurate 
knowledge of the stellar radii of the host stars. 

This problem is particularly acute for the planetary 
system around KOI-961. Because the host star is an 
M dwarf, the photometry-based stellar properties in the 
Kepler Input Catalog suffer from large systematic errors. 
However, we have taken advantage of two opportuni- 
ties afforded by the nature of this star. First, we find 
that KOI-961 is a near twin to a nearby, well-studied M- 
dwarf. By performing a differential analysis with respect 
to Barnards star, we find that KOI-961 is also a small 
star with large proper motion. Second, the high proper 
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Fig. 14. — The probability distribution for the radius of KOI 961.01, accounting for potential binary star configurations of KOI 961 that 
could introduce diluting flux within the Kepler aperture. Scenario (a) is that KOI 961 is a single star (probability 0.76), scenario (b) is that 
KOI 961 is a binary with the planet transiting the brighter component (probability 0.12), and scenario (c) is that KOI 961 is a binary with 
the planet transiting the fainter component (probability 0.12). The planets 961.02 and 961.03 have analogous probability distributions. 

and transit depths for all three planets. The transit 
depths combined with our estimate of the stellar radius 
yield planet radii that are all less than the radius of the 
Earth. We measure radii of 0.78 ± 0.22 R ffi for KOI 
961.01, 0.73 ± 0.20 R® for KOI 961.02, and 0.57 ± 0.18 
R ffi for KOI 961.03. 

Incorporating all available observational constraints on 
the presence of potential blends (both chance- alignment 
and hierarchical), we are able to establish that the tran- 
sit signals are indeed planets in the KOI 961 system. 
We also quantify the probability that KOI 961 might be 
a binary with planets transiting either of the two com- 
ponents (which would lead to an underestimate of the 
planet radii) , and find that this probability is low enough 
to remain confident that we have not drastically misclas- 
sified the planets. Follow-up adaptive optics imaging of 
KOI 961 will be able to put even tighter constraints on 
this possibility. 

As a planet-hosting star, KOI 961 provides an impor- 
tant data point for the study of the planet-metallicity 
relationship in the context of low-mass stars and 
terrestrial-sized planets. The planet-metallicity cor- 
relation was originally discovered for Sun-like F, G 
and K-type stars, and correlated the stellar metal- 
licity to t he likelihood of hosting a Jovian-mass 
planet (e.g. iSantos et all 120011: iFischer fc Valentil l2005t 
IJohnson et al.ll2010al) . Recently developed methods for 
estimating the metallicities of M dwarfs have shown 
that M-dwarf planet-hosts are also systematically metal- 
rich (IJohnson fc Adds! 120091: iRoias-Avala et al.1 120101 : 
iSchlaufman fc Laughlinl 120101 : iNeves et al.ll201lD . Being 
significantly metal-poor, KOI 961 represents an outlier 
in the context of these studies, but it is not clear that 
the planet-metallicity correlation extends to stars which 
host Earth-mass and sma ller exoplanets. Recent analy- 
sis by lMavor et al.1 ()2011[ ) of radial- velocity detected ex- 
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Fig. 15. — Planetary Mass vs. Radius showing the ±1<t values 
for the KOI 961 planetary radii and the range of possible masses. 
We show theoretical models for cold spheres composed of pure 
iron, rock (Mg2SiQ4 ) and water-ice, interpolated from Table 1 of 
IFortnev et al. (20Q7J). We also include Solar System bodies Mer- 
cury, Mars, Venus and Earth, from left to right (black triangles). 
Planetary bodies are unlikely to contain significant amounts of el- 
ements that are heavier than iron, so we use the pure iron mass- 
radius contour to place upper limits on the masses KOI 961.01 
(red), .02 (green) and 0.3 (blue). 



motion greatly facilitates the validation of the planetary 
system based on examination of historic imaging data. 
We find that there are no potential eclipsing binaries at 
the current position of KOI-961 up to 4.8 magnitudes 
fainter. 

Using revised and refined stellar parameter for KOI- 
961 we perform a reanalysis of the transit light curves, 
confirming the three signals at the periods reported by 
iBorucki et al.l f|201 lh . We revise the impact parameters 
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oplanets indicates that the correlation may not apply to 
stars which only host planets less-massive than Neptune. 
However, there is evidence that the planet-metallicity 
correlation is important for small-planets around low- 
mass stars, based on a comparison of the KOIs to Ke- 
pler field stars by iSchlaufman fc Laughlinl (1201 lb. The 
discov ery of KOI 961, and the studies o f iMavor et al.l 
(|2011h and ISchlaufman fc Laughlinl (|2011[ ). strongly mo- 
tivate further investigations of the planet-metallicity cor- 
relation for low-mass planets around low-mass stars. 

The detection of the planetary system around KOI 
961 is remarkable since it is one of only a few dozen 
mid-to-late M dwarfs currently observed by Kepler. Fig- 
ure [16] shows the r-K distribution for a ll targets cur- 
rently observed by Kepler, and KOIs in iBorucki et al.l 
(|2011h . We have removed all objects with J - K > 
1.0, as these are likely distant giant stars with signif- 
icantly reddened colors (se e Figure HI). KOI 96 1 is the 
reddest dwarf KOI from the Bo rucki et al.l 1)20111 ) release, 
and is also likely the least massive KOI. Only 27 dwarf 
stars observed by Kepler have r-K greater than KOI 
961 's value of 4.444. Combined with the low probability 
of a planetary system being geometrically aligned such 
that a transit is observed (only 13% in the case of KOI 
961.01), Kepler's discovery of planets around KOI 961 
could be an indication that planets are common around 
mid-to-late M dwarfs, or at le ast not rare. This wo uld be 
consistent with the results of iHoward et al.l (|2011l ). who 
used the Kepler detections to show that the frequency 
of sub-Neptune-size planets (Rp = 2-4 R ffi ) increases 
with decreasing stellar effective temperature for stars ear- 
lier than MO. Results from on-goin g exoplanet survey s 
of M dwarfs such as MEarth (e.g. Ilrwin et all 1201 lbl ). 
and future programs such as the Habitable Zone Planet 
Finder (e.g. iMahadevan et al.ll2010HRamseY et al.ll2008H 
and CARMENES (jQuirrenbach et all l2010f h will shed 
light on the statistics of low-mass planets around mid- 
and later M dwarfs. 
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